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h^- , An exactly solvable variant of mixed spin-(l/2,l) Ising-Heisenberg diamond chain is considered. 

Vertical spin-1 dimers are taken as quantum ones with Heisenberg bilinear and biquadratic interac- 
tions and with single-ion anisotropy, while all interactions between spin-1 and spin-1/2 residing on 

£T^ , the intermediate sites are taken in the Ising form. The detailed analysis of the T = ground state 

phase diagram is presented. The phase diagrams have shown to be rather rich, demonstrating large 
variety of ground states: saturated one, three ferrimagnetic with magnetization equal to 3/5 and 
another four ferrimagnetic ground states with magnetization equal to 1/5. There are also two frus- 
trated macroscopically degenerated ground states which could exist at zero magnetic filed. Solving 
a the model exactly within classical transfer-matrix formalism we obtain an exact expressions for all 

thermodynamic function of the system. The thermodynamic properties of the model have been de- 
scribed exactly by exact calculation of partition function within the direct classical transfer-matrix 
formalism, the entries of transfer matrix, in their turn, contain the information about quantum 
states of vertical spin-1 XXZ dimer (eigenvalues of local hamiltonian for vertical link). 
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I. INTRODUCTION 



e 

Lattice models of quantum magnetism continue to be in the focus of attention of theoretical condensed matter 
physicists. Beside of great practical importance connected with description of magnetic and thermodynamic properties 
of real magnetic materials, this research area is also attractive from the general statistical mechanics and strongly 
correlated system theory points of view, especially when one deals with exactly solvable strongly interacting many-body 
system. The diamond chain is a one-dimensional lattice spin system in which the vertical spin dimers alternate with 
single spins(See figure [T]). This model with S — 1/2 is believed to describe the magnetic lattice of mineral azurite, 
CU3 (003)2(011)2, which is famous for its deep blue pigmentation[j}jj| . Theoretical research of various aspects of 

0^ ■ diamond chain physics received much attention during last years[5|-[l6|. Diamond chain and especially diamond chain 

with mixed (S, S/2) spin are shown to have very rich ground state phase diagrams with Haldane and several spin- 

^— . , cluster states which are tensor product of exact local eigenstates of cluster spins [a, @. Many other issues of diamond 

chain physics have been investigated theoretically during last years including Dzyaloshinsky-Moriya term influence on 
magnetization processes[7], multiple-spin-exchange effects [8], magnetization plateaus [9J,[lO|, magnetocaloric effect |ll|. 
e.t.c. Very recently another one interesting feature of diamond chain, possibility of localized magnon excitations, has 
been also investigated [12]. 

Especially important issue is the effect of frustration, which is rather strong in antiferromagnetic diamond chain 
because of triangular arrangement of the sites. Variants of frustrated recurrent lattices with diamond plaquette have 
been also studied in Refs. [21] and [22]. However, diamond chain is not integrable in general. Thus, the exact analysis 
of the dynamic and especially thermodynamic properties of diamond chain is a very complicated issue. Nevertheless, 
one can consider various exactly solvable variants of spin systems possessing diamond chain topology of interaction 
bonds with simplified structure of interactions [13|- [16]. The diamond chain with only Ising type of interaction have 
been exactly solved within classical transfer-matrix technique in Ref. [l3| revealing reach T = ground state phase 
diagram. Yet another exactly solvable diamond chains have been considered in Refs. (|14]-[17|). where vertical spin 
dimers have been taken as quantum ones with XXZ interaction while interaction between spins localized on vertical 
dimer sites and spins from the single sites alternating with them is of Ising type. Mapping the system into a single 
Ising chain within iteration-decoration transformation |l8l420l | the authors gave complete description of the ground 
state properties, T — ground state phase diagram as well as thermodynamic functions for all values of vertical 
dimer spins magnitude S. In a very recent paper the diamond chain with XX-interaction has been considered in the 
Jordan- Wigner formalism |17|. 
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Figure 1: The diamond chain with Ising and Heisenberg bonds, 
correspond to Ising interactions. Large(small) circles denote S (a) 



Solid bold lines denote XXZ quantum bonds, dotted lines 
spins. 



Considering the mixed spin chains (or another one-dimensional spin systems with more complicated geometry) with 
Ising and Heisenberg bonds (or even just Ising counterparts of the known quantum spin models) one can achieve two- 
fold goal: to construct novel exactly solvable lattice spin model which allow one to obtain analytic expression for all 
thermodynamic functions of the model, and to get a n appro ximate models which can be useful for understanding the 
properties of underlying purely quantum models [14H161 |23j L 29] . Exact thermodynamic solutions of such models even 
can shed a light to the properties of real magnetic materials. For instance, for alternating spin chain even simplest 
models with only Ising interaction can reflect the underlying magnetic behavior of the corresponding Heisenberg 
counterpart at least in qualitative way [30l . l3ll | , moreover, some exactly solvable models with Ising and Heisenberg 
bonds can also provide satisfactory quantitative picture[29j. Very recently, the synthesis of novel class of trimetalic 
3d-4d-4f coordination polymers has been reported. One of them, a Id coordination polymer compound containing 3d 
(Cu 2+ ), 4d (Mo 5+ ) and 4f (Dy 3+ ) ions is shown to exhibit the properties of Id magnet with Ising and Heisenberg 
bonds [32]. The appearance of Ising interactions between magnetic ion in this compound is connected with the extremely 
anisotropic properties of Dy 3+ ground states (g\\ = 19.6, g± « 0). Thus, the interactions of Dy 3+ with the surrounding 
Cu 2+ and Mo 5+ ions are, to the great extent, of Ising type, involving spin projection along the dysprosium anisotropy 
axis. While, the interaction bonds Cu 2+ — Mo 5+ correspond to Heisenberg interaction[32j. Though, the aforementioned 
Id coordination polymer system is not exactly the Ising-Heisenberg diamond chain considered in the present paper, 
this discovery of novel classes of magnetic materials makes the investigation of exact solutions of spin chains with Ising 
and Heisenberg bonds important form practical point of view as well. 

In this paper we consider mixed spin- (1,1/2) diamond chain with Ising and Heisenberg bonds which extends the 
system considered in Ref. [16j by including biquadratic term for S = 1 XXZ-dimers and single-ion anisotropy. 
Biquadratic terms are usually originated from the spin-lattice coupling in the adiabatic phonons approximation[33| 
but also can be considered as the effect of quadrupole interaction between the spins. We do not make any assumption 
about the origin of biquadratic terms considering them as a part of general Hamiltonian. The model allow one 
to calculate the partition function and, thus, all thermodynamic quantities exactly within classical transfer-matrix 
formalism |34j. We present the analysis of T = ground-state phase diagrams and plot the curves of magnetization 
processes for finite temperatures, demonstrating magnetization plateaus at M = 1/5 and M = 3/5 in the units of 
saturation magnetization. 

The paper is organized as follows. In the second Section we formulate the model and present eigenvalues of block 
Hamiltonian. In the third section we describe possible ground states of the system and present the ground-states phase 
diagram. In section 4 we present its exact solution and discuss the magnetization and thermodynamics properties of 
the model. Finally in sec. 5 a short summary is followed. 



II. THE MODEL AND ITS EXACT SOLUTION 

Let us consider the system of vertical S — 1 spin dimers with Heisenberg XXZ bilinear and biquadratic interactions 
and uniaxial single-ion anisotropy. These dimers are assembled to the chain by alternating with Ising spins <j, so that 
each spin S in certain dimer interacts to its both left and right Ising spins via Ising type interaction (See figure [T]). 
So, we have the so-called diamond-chain with 5 = 1 Heisenberg dimers and a = 1/2 Ising spins between them. The 
corresponding Hamiltonian is the sum over the block Hamiltonians: 
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U h = J(8u ■ S l2 ) A + K(S n ■ S i2 )i + D ((S^) 2 + (S z l2 ) 2 ) + Jbfa + <T i+ i){Sl x + S z l2 ) - h x {Sl x + S£), 



where N is the number of the unit cells (blocks with two spin-1 and one spin-1/2), while i correspond the particles at 
i-cell. J being the coupling constant of bilinear AAZ-interaction between Heisenberg spins, which we assume to be 
of the following form 

(S;i • S; 2 )a = A(Sf l Sf 2 + Sn^v) + Si X Sf 2 , (2) 

whereas K denotes the biquadratic XXZ-interaction term, D means the single ion-anisotropy and Jo being the purely 
Ising interaction term, hi and h\ are the external magnetic field acting on cr, and Si respectively. 

In order to solve this model, at first, we need to diagonalize the block Hamiltonian for arbitrary i-th block. Nine 
eigenvalues of %$, A n (cri, o"i+i), n — 1...9 can be found analytically, which write down as 

Ai, 2 = J + K + 2D±2[-hi + J (<ri + a i+ i)], 

X 3A = A(J + AK)+D±[-h 1 + J (a i + a i+1 )}, 

A 5i6 = -A{J-AK) + D±[-h 1 + J (a l +a l+1 )], 

A 7 = -J + K + 2D, 

-J+{l + AA 2 )K + 2D 1 
A s,9 - ^ ± 2 ' ^ ' 

where for simplicity R denotes the following expression,, 



R = y/8A 2 (J - K) 2 + (J-K- 2D) 2 . (4) 

The eigenvectors of block Hamiltonian Hi up to the inversion of all spins are 

l«a>=|M>. => 
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where 1,0 and —1 in first(second) place stand for the S z = 1,0 and —1 states for first(second) spin in vertical dimer 
and the following notation is adopted: 

1 / 2D±R\ 

c± = 2a{ 1 + 1<-t)- ^ 

The eigenvectors of the dimer defined by \vi) corresponds to the parallel ordered spins with magnetization per 
site m s = 1, the eigenvectors |«i, s ) and |wi jQ ) correspond to symmetric and antisymmetric states vector respectively, 
whereas \vo, a ) is an antisymmetric state vector with magnetization m s = 0, finally |uo,±) are the symmetric eigenvectors 
with magnetization m s = 0. 

The remaining eigenvectors of the eigenvalues A2, A4 and Xq can be obtaining using the spin inversion. 

A. Special case K = J 

At the special value K = J the S* ot = sector of the block Hamiltonian undergoes qualitative changes which can 
be obtained substituting carefully K = J value into the general solution presented above in eq.((S]). This should be 
considered as a consequence of the special symmetry of the Hamiltonian for these values of parameters, more precisely, 



the fact that operator (Si • S2)a + (Si • 82)^ can be represented in terms of permutation operators P12. In this case 
the eigenstates |i>o,a) and \vo,±) 0I the Hamiltonian reduce to the following ones 

|«d,a>=^(|l,-l>-|-l,l», => A 7 = 2A 

K+)=|0,0), => A 8 = 2JA 2 , 

K_)=^=(|l,-l) + |-l,l)), => A 9 = 2(JA 2 + £). (7) 

Note that a straightforward substitution in eq.© could yield to an undefined coefficients of the eigenstates. 

III. GROUND STATES PHASE DIAGRAMS 

Let us describe possible T = ground states of the chain under consideration and the corresponding energies per 
block. Generally speaking, there are 9 x 2 = 18 possible ground states for each block. However, if one restricts 
himself with the ground states which are equivalent up to the inversion of all spins one will arrive at the following spin 
configurations. The fully polarized state (M = 1, here M is magnetization per spin) 

N 

\SP)=\{ |«2>i <8> I t)i, e S P = J + K + 2D + 2J Q -%H, (8) 

t=i 

where | f)i(| i)i) stands for the up(down) state of the cr-spin in the i-th block. Hereafter, we also put hi=h,2=H ■ The 
next sector of ground states contains three different ferrimagnetic spin configurations with the value of magnetization 
equal to 3/5 (M = 3/5), 

N 

\Fl) = I] M* ® I i)h e F1 = J + K + 2D- 2J - -H, 

i=\ 

N 

\F2) = J] \v ltS )i <g> I t>t, £F2 - A(J + AK) + D + J - -H, 

1 = 1 

N 

\F3) = n K«>< ® I t)i, e F3 = -A(J - AK) + D + J - -H. (9) 

i=i l 

There are also four another ferrimagnetic ground states with M — 1/5: 

N 1 

1^4) = J] \v ltS )i ® I l) h e F4 = A(J + AK) + D- J - -H, 

N 1 

1^5) = J[ \v 1>a )i ® I i) i; e F5 = -A(J -AK) + D-J Q - -H, 

i=l 

W 1 

1^6) = JJ \vo, a )i ® I t)i, £F6 = -J + K + 2D- -H, 

i=i 2 

W 1 1 

1^7) = JJ K-)i ® 1 1)<, £f 7 = ^(-^+(l + 4A 2 )if + 2 J D-i?)--ff. (10) 
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When no external magnetic field is applied, there are also possibility of frustrated ground state formation, in which 
the orientation of a spins in each block are not defined. In this case all a spins become decoupled and behave like free 



spins. There are two frustrated ground states 

N 



|FF1) = J] | VOj _), ® |0i, £^2 = ~ (- J + (1 + 4A 2 )F + 2D-R), 

i=\ 

N 

\FR2) = J] \v 0>a )i <g> |0<, era = - J + A 1 + 2£>. (11) 



Here |£).; stands for arbitrary value of the a spin in the i-th block. Thus, for H = if S = 1 dimer has Sf ot = then 
its neighboring er-spins become decoupled which results in macroscopic non-zero entropy S/N = log 2 for each of the 
frustrated ground states of Eq. 1111 Applying magnetic field one removes the two-fold macroscopic degeneracy driving 
\FR1) and \FR2) ground states into |F7) and |F6) respectively. Nevertheless, in some papers two last non-degenerated 
ground states of eqs. fTU)) are mentioned as frustrated ones [Ij-Qjj]. 

Hereafter, to discuss the phase diagrams we will consider the external magnetic field as hi=fi2=H. It is also 
convenient to present all parameters in the units of \J\. Thus, we define k = K/\J\, jo = Jo/\J\, S = D/\J\ and 
h = H/\J\. In Figure [2] one can see four ground states phase diagrams plotted in (K,/i)-plane demonstrating vast 
variety of ground states for fixed values of the 5 and jo- These plots summarize the effect of biquadratic term. The 
equations of phase boundaries for J > 0, A = 2, jo — 0.5 and 5 = 0.5(Figure[5]Ja)) are 

between F7 and SP, h =- f 6 - 15k + ^k 2 + 32(k - 1)A , (12) 

between F7 and F3, h =- (-2 - 9k + ^ n 2 + 32(k - l) 2 

between F3 and SP, h =4 — 3k. 

The rest phase boundaries in this case are either horizontal or vertical lines in (k, /i)-plane. F5 and F3 as well as -F6 
and Fl are separated by the h = 1 line, while the phase boundary between F5 and F6 as well as between F3 and F\ 
are vertical line is situated at k = 2/3. Ground states Fl and F5 are separated by the line k = 1/12(— 17 + V337) ~ 
0.11313. 

For the case of antiferromagnetic Heisenberg interaction between 5 = 1 spin (J < 0) presented in Figure [5{b) the 
equations of the phase boundaries are 



between F7 and SP, 
between F7 and F2, 



=- (-15k + v/32(k+1) 2 + (k + 2) 2 ) , (13) 

■■- (-4 - 9k + v/32(k + 1) 2 + (k + 2) 2 ) , 



between F2 and SP, h =2 - 3k, 

between F4 and Fl, h =1 — 3k. 

The boundary between Fl and saturated ground state SP is the straight line h = 2. Horizontal line k = separates 
FA and Fl as well as Fl and F2. In order to demonstrate the significant role of exchange anisotropy A we plotted 
also the ground states phase diagrams for the isotropic case A = 1 ( Figure [^c) and (d)) for antiferromagnetic and 
ferromagnetic J respectively. One can see simplification of the ground state phase diagram via the disappearance of 
2 ground states presented in the case anisotropic case A = 2( Figure^a) and (b)). So, for J > one can find, beside 
SP, only F3, F5, F7 and for J < only Fl and F7 ground states respectively. The equation of phase boundaries for 
antiferromagnetic J, isotropic A = 1 and jo — 0.5, 6 — 0.5 are: 



between F7 and SP, h 



i(6-3K+V8(K-l) 2 + K 2 ), (14) 

i(-K +v /8(K-l) 2 + K 2 ), 



between F7 and F3, h 

2 

Ground states F3 and SP are separated by the horizontal line h — 3, another two straight lines appears between F3 
and F5 and F7 and F5 at h = 1 and k = 1/7 respectively. In the case of ferromagnetic J and isotropic exchange 
interaction presented in Figure HJd) one can see only three possible ground states Fl, F7 and SP. Here the horizontal 
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Figure 2: Ground states phase diagrams in (K,ft)-plane demonstrating the effect of biquadratic term. Here k = K/\J\ and 
h — H/\ J\ The values of all other parameters, jo = Jo/\J\, 5 — D/\J\ are fixed as follows: jo = 1/2, S — 1/2. Left panels ((a) 
and (c)) corresponds to the antiferromagnetic Heisenberg interaction J > 0, the right panels ((a) and (c)) - to ferromagnetic 
one J < 0. In upper panels ((a) and (b)) A = 2 has been taken, while in lower panels ((c) and (d)) one can see phase diagrams 
for isotropic case A = 1. 



line h — 2 separates Fl and SP while other two phase boundaries are given by 

between F7 and SP, h=-U-2,n+ y/8(n + l) 2 + re 2 ) , 



between F7 and Fl, 



(15) 



h =- (2 - 3k + v / 8(k+1) 2 + k 2 ) , 



In order ro summarize the effects of the Ising coupling Jo we plotted another two ground state phase diagrams 
presented in Figure [3] Here, the left panel shows ground states boundaries in (jo, /i)-plane for fixed values of 5, k and 
A, while right panel demonstrate the phase boundaries for fixed S, h and A in (re,jo)-plane. For the sake of briefness 



we just list the equation of phase boundaries. For left panel: 

between F3 and SP, h =-3- + jo, (16) 

8 

between F3 and F7, ft=-( — ] + jo, 

between F3 and F5, 
between F5 and Fl, 
between Fl and SP, 

between F7 and F5, 

between F3 and Fl, 
For right panel: 

between F7 and SP, j = - - (l + ^2{k - l) 2 + (n + 1)A , (17) 

between F7 and F3, j =- ( 2- -« - ^2{n - l) 2 + (k + l) 2 J , 

between F7 and F5, j =- ( --« + \J2(k - l) 2 + (n + l) 2 J , 
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between F7 and Fl, j =- f 3 + ^(k - l) 2 + (k + l) 2 

3 / 1 
between SP and F3, jo = I 1 H — K 

between F3 and F5, jo =-, 

between F5 and Fl, jo =— I 1 H — ft 



2 V 2 

Another two ground state phase diagrams demonstrating the influence of single-ion anisotropy are presented in Figure 
|U Left(right) panel exhibits phase boundaries for for h = 1, jo = 1 and A = 0.5 (k — 0.5) respectively. The phase 
boundaries for the left and right panels respectively are 

3(2 + 2k — n 2 ) 

between Fl and F7, 5 = - -—-: r—^, (18) 

4(2 + K) 

between Fl and F5, S =-(2 — «), 



4 



between F5 and F7, 



6 =- (l - « + v / 2 + 10k+1/4:kA 



and 



between Fl and F4, 5=-(A + l) 2 , (19) 

between Fl and F5, 5=-(A-l) 2 , 

between F7 and F4, 6 =- fl/2 - ^9/4 + A 4 - 4A 3 + 5A 2 - 6A 

1 



between F7 and F5, 5 =- ( 1/2 + ^9/4 + A 4 + 4A 3 + 5A 2 + 6A 



(20) 
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Figure 3: Ground states phase diagrams demonstrating the influence of Jo. Here A = 1/2, k = 1/2 and 5=1. The panel (a) 
shows the ground states phase diagram in (jo,/t)-plane for fixed value of biquadratic interaction k = 1/2; the panel (b) shows 
the ground states phase diagram in («, jo)-plane for fixed value of magnetic field h — 1. 

And finally, the zero field ground states phase diagram are presented in Figure [5j The phase diagram exhibits the 
appearance of two frustrated ground states at zero magnetic field. We choose D — and rather small value of 
biquadratic interactions K = 0.1. For convenience we normalize parameters by Jo > 0. The equations for the phase 
boundaries are 



between Fl and F5, 
between Fl and FR2, 

between FR2 and F5, 
between F5 and FR1, 

here J = J/Jq and K = K/Jq. 



J = 



1-K{1-A 2 ) 



(21) 



J=l, 



J = 



1 + K(l- A 2 ) 



2A - 1 + A(l + A)(l + 2A)JsT + (1 + A(l + A)K)Vl + 8A 2 
~~ 2A(1 + A) 



(22) 



IV. EXACT SOLUTION AND THERMODYNAMICS 



The present model could be solved exactly using the known decoration transformation early presented by M. E. 
Fisher [la. [l9J. and recently generalized for arbitrary spin interactions [20], where one maps the partition function of 
the system to the partition function of Id Ising model, writing down the relations for the entries of transfer matrix 
and thus obtaining the relations between model parameters and that of Ising chain. But here we implement a direct 



1 



-1 

-2 
-3 
-4 
-5 



I ■ I ■ LJ— 4 ■ I *=1 5 




(a) 



i . i . i . i . i 



■10 12 3 4 
K/J 




Figure 4: Ground states phase diagrams demonstrating the effect of single- ion anisotropy. Here h — 1, jo = 1. Panel (a) exhibit 
the ground states phase diagram in (k, 5)-plane for fixed value of exchange anisotropy A = 1/2; panel (b) exhibits the ground 
states phase diagram in (A,5)-plane for fixed value of biquadratic interaction k = 1/2. 

transfer matrix calculations without any account to the solution of Ising chain. Therefore let us consider the following 
partition function of the system, which can be represented as 



N 



2 



Z = ^2 S Ps exp(-/3%) =^ JJexp(/3 — {a l + <r l+1 ))Z(a t , <r i+1 ), 

a a i—1 

where f3 as usual is inverse temperature. Here 



(23) 



Z(a il a i+ i) = Sp i cxp(-/3H l ) = J^ exp(-/?A„(cr l , cr l+ i)). 

n=l 
Then, using ([3]) one can express the one-block partial partition function Z{a i ,<r i+ i) in the following form: 

2 

Z(ai,<T i+1 ) = 22 Z n cosh (/3n(hi - J (ai + a i+ i))), 

n=0 

Z a = e^ J - K - 2D ^ + 2e^ J ^ WA ^ K - 2 ^cosh(^), 

Z x = 4e-' 3 ( AK + I5 'cosh(/3AJ), 

z = 2e -p(J + K+2D)_ 



(24) 



(25) 



After that, the partition function f|23[) take the form similar to the partition function of a chain with classical two 
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Figure 5: Zero field ground states phase diagram drawn in the (A, J/Jo)-plane for fixed D = and small values of K/Jo = 0.1. 
The frustrated ground states are exhibited. 



state variables on each site: 



N 



where Ai 2 are two eigenvalues of the transfer-matrix 

_ ( e^Z+ Z \ 



where 



Z± = Z (±1/2, ±1/2), 

Z = Z(l/2, -1/2) = Z(-l/2, 1/2). 



(26) 



(27) 



(28) 



Then, calculating the eigenvalues and taking thermodynamic limit, when only the largest eigenvalue survives, one 
arrives at the following expression for the free energy per block: 



/ = -ilogQ(, ^5-+ t -^5 +yU^2_. -, ^5_)^+4Z„-- 



(29) 



As soon as the free energy per block is calculated, one can obtain analytic expressions for all thermodynamic function. 



A. Magnetization and quadrupole moment 



Magnetic quantities can be obtained using the free energy expression obtained in (1291) . Therefore the magnetization 
of the spin-S* can be expressed as 
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Figure 6: (a) Magnetization as a function of K/J, for H/J — 3.0, T/J — 0.15, A = 2 and Jo/ J — 0.5. (b) Magnetization 
as a function of H/J, for if/ J = 0.5, T/J — 0.15, A = 1 and Jo/ J = 0.5. (c) Magnetization versus temperature T/J, for 
Jo/ J = 0.5, H/J = 4, A = 1.4 and D/J = 1 . 



M s = 



27VZ 



H S Ps 



while the magnetization of spin-er reads 



1S& + Sf 2 )e 
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Thus, the total magnetization is given by 
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(32) 



Figure Ela) displays the plot of magnetization as a function of biquadratic interaction K in units of J , for fixed 
values of H/J = 3, T/J = 0.15, A = 2 and Jq/J = 0.5. Here one can see the plateaus for different values of J, these 
plateaus occur as expected at 1/5 and 3/5. The plots of magnetization processes (M versus H) for the system under 
consideration is presented in figure E^b). Here the values of parameters are fixed as K/J — 0.5, T/J = 0.15, A = 1 
and the value of D/J varies from -1 to 1. As in can see from the phase diagram in Figure ©(c), two intermediate 
magnetization plateaus at M = 1/5 and M = 3/5 exhibited by the system in the case correspond to F3 and -F5 
ground sates for T — 0. Thermal behavior of magnetization at the fixed external field is presented in figure [(^c) where 
the following values of parameters are assumed: H/J = 4, A = 1.4, D/Jq = 1 and Jq/J — 0.5, the competition 
between ferromagnetic state SP and ferrimagnetic state F\ with magnetization M = 3/5 is occurred with increase of 
the temperature. Close to zero temperature one obtains three well defined values for the magnetization, which are in 
accordance with plateaus displayed in figure [BJa-b). 

As the system under consideration contains sites with spin-1 one can define another important physical quantity, 
quadrupole moment, which can be obtained by the thermodynamic relations as well 



Q = 



2NZ ^ 



Sp s 



;((^) 2 + (5f 2 ) 2 )e-^ 



df_ 

dD 



(33) 



PMM 



In figure Eta) the plots of quadrupole moment as a function of K/J for fixed values of H/J = 3, Jq/J = 0.5, 
A = 2 and T/J = 0.15 are presented for several temperatures. The non-trivial and non-monotone behavior of Q under 



12 



^r-O.l — = 0.15 y 



K -1 


K 


1 


A" 


1 


•*" -1 


J ' 


J 


2 


J 


2 


./ ' 





™V\ 








\ \ 








\ \ 






(c) 


\ s 


•^ 












^^"^-^^^I 








-— ~* - 











Figure 7: Quadrupole Moment: (a) As a function of K/J, for .ff/J = 3, Jo/ J = 0.5, A = 2 and T/J = 0.15. (b) As a function 
of D, for J/ Jo = -0.5, T/J = 0.1, 7J/J = 2 and A = 2. (c) As a function of T, for J /J = 0.5, D/J = 1.0, F/J = 4 and 
A = 1.4. 

variation of K can be understood if one take into account appearance of F6 and FR2 ground states in which vertical 
quantum dimer is in |i>o,— ) eigenstate. Calculating expectation value for the operator Q for this state one obtains: 



1 



(«b,-b((5i*r+(5inK_) 



i 



i + ^- (i i 2d - r \ 

1 + 8A 2 ^ + K-J ) 



(34) 



which actually defines low-temperature behavior of the quadrupole moment. The quadrupole moment dependence of 
the uniaxial single-ion anisotropy parameter D is illustrated in figure[7Ih) for H/J — 2, Jo/ J = —0.5, T/J = 0.1 and 
A = 2. Similar to the magnetization case, we obtain some plateaus, but for higher values of D we have a decreasing 
curve instead of plateaus(solid line). On the other hand, as soon as the temperature increases these plateaus obviously 
disappear. In the Figure Etc) we show the quadrupole moment as a function of the temperature for fixed values of 
Jq/J — 0.5, D/J — 1.0, H/J = 4 and A = 1.4, similar to the case of magnetization (figJBJc)), the quadrupole moment 
leads to fixed values at low temperature which are related to the plateaus displayed in figures [TJa-b) , while at high 
temperature average quadrupole moment leads to 2/3, which corresponds to equal probabilities for all three values of 
5=1 spin. 

B. Entropy and specific heat 

The entropy of the system can be obtained according to general thermodynamic relation, 



S(T) = - 



d£ 
dT 



(35) 



hi,h z ,D 



One can see the plot of entropy S(T) of the model as a function of K/J in figure [SJ here we fix the following values 
of the parameters: Jo/ J = 0.5, T/J = 0.1, A = 3 and D/J = 1, assuming the external magnetic field is fixed at 
H/J = 1.0, 0.5 and 0.1. Three panels correspond to three values of A, (a) A = 0.5, (b) A = 1, (c) A = 3. A series of 
peaks corresponding to quantum phase transition points appear on the curve which is in accordance with the phase 
transition at zero temperature illustrated in figure [5] According to previous discussion, when external magnetic field 
vanishes there are a residual entropy of the model corresponding to frustrated ground states. This residual entropy 
is a constant value 5(0) = ln(2)=0.693 related to the ground state degeneracy. For particular values of the magnetic 
field like H/J — 1 it is possible to occur a twofold degeneracy of the ground state energy. 

Using the entropy S(T), one can obtain the specific heat 



<-<§ 



(36) 



H 



13 



H 




H 




H 






= 1.0 — 




= 0.5 -- 




= n 


J 




J 




J 






Figure 8: The entropy SiT) against J: assuming fixed values for J /J = 0.5, T/J = 0.1, D/J = 1. (a) A = 0.5, (b) A = 1, (c) 
A = 3. 
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Figure 9: The specific heat C(T), for the fixed values of J /J = 0.5, T/J = 0.2, J /J = -0.5 and D/J = 1. (a) A = 0.5, (b) 
A = l, (c) A = 3. 

In order to study the specific heat properties for the model under consideration, in fig[9] we display the specific heat 
as a function of the parameter K/J for the fixed values of Jo/ J = 0.5, T/J = 0.2 and D/J = 1. For three values of 
A, (a) A = 0.5, (b) A = 1, (c) A = 3. In figure |H] the corresponding plots exhibit four peaks which are related to 
the phase transition at zero temperature. As one can see the phase transition effects appears for small values of the 
biquadratic term K . 



CONCLUSIONS 



We have considered an exactly solvable variant of diamond chain with mixed S = 1 and S = 1/2 spins. The 
vertical S = 1 dimers are taken as quantum ones with Heisenberg bilinear and biquadratic interactions and with 
single-ion anisotropy terms, while all interactions between S = 1 spins and S — 1/2 spins residing on the intermediate 
sites are taken in the Ising form. The model generalizes the model of diamond chain with Ising and Heisenberg 
bonds considered in (29|. Our results supplement the previously obtained ones for the case of S = 1 vertical XXZ- 
dimers with only bilinear Heisenberg interaction. The detailed analysis of the T = ground state phase diagrams is 
presented as well as the exact plots of various thermodynamic functions. The effect of biquadratic term and single- 
ion anisotropy are summarized in the corresponding phase diagrams. The phase diagrams have shown to be rather 
rich, demonstrating large variety of ground states: saturated one, three ferrimagnetic with magnetization equal to 
3/5 and another four ferrimagnetic ground states with magnetization equal to 1/5. There are also two frustrated 
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macroscopically degenerated ground states which could exist at zero magnetic filed. The thermodynamic properties 
of the model have been described exactly by exact calculation of partition function within the direct classical transfer- 
matrix formalism, the entries of transfer matrix, in their turn, contain the information about quantum states of vertical 
S = 1 XXZ dimer (eigenvalues of local hamiltonian for vertical link) . The plots of entropy and specific heat are also 
presented. 
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